Excitation of coherent polaritons in a two-dimensional atomic lattice 
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We describe a new type of spatially periodic structure (lattice models): a polaritonic crystal 
(PolC) formed by a two-dimensional lattice of trapped two-level atoms interacting with quantised 
electromagnetic field in a cavity (or in a one-dimensional array of tunnelling-coupled microcavities) , 
which allows polaritons to be fully localised. Using a one-dimensional polaritonic crystal as an 
example, we analyse conditions for quantum degeneracy of a low-branch polariton gas and those for 
quantum optical information recording and storage. 
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q ■ I- INTRODUCTION 

O: 

This work develops the ideas that S.A. Akhmanov paid much attention to in the last years of his life. It addresses 
the general principles, laid by Akhmanov, behind the optical information recording and processing using non-linear 

i 1 imaging in spatially periodic or inhomogeneous dynamic structures excited by laser radiation in nonlinear media of 

i-Ch ' various kinds (see, e.g., [1]). 

In recent years, great advances have been made in laser control of macroscopic amounts of ultracold atoms [2]. The 
ability to produce an array of macroscopic atomic Bose-Einstein condensates (BECs) by cooling and trapping atoms 
in one- and two-dimensional optical lattices enables research into various physical aspects of phase transitions. Strong 
atom-photon coupling has recently been demonstrated for BEC atoms in a cavity [3]. Basically, the fabrication of 
q-i spatially periodic atomic structures confined in an optical cavity opens up novel opportunities for gaining insight 
into critical phenomena in coupled atom-photon systems [4]. Recent advances in nanotechnology and photonics 
have made it possible to create such structures using one-dimensional (ID) arrays of coupled microcavities, so-called 
coupled- resonator optical waveguides, containing two- or three-level atoms [5-7]. A key role in determining the 
behaviour of such systems is played by bright- and dark-state polaritons-bosonic quasi-particles resulting from a 
linear superposition of a photon and macroscopic (coherent) excitation of a two-level atomic system. 
^sO ■ In Ref. [8] Bose-Einstein condensation and the Kosterlitz-Thouless phase transition for polaritons resulting from 
the interaction of a quantised light field with an ensemble of two-level atoms in a cavity is proposed. Note that 
the phase transition in question may take place at sufficiently high (room) temperatures because of the low-branch 
polariton effective mass [9]. Kasprzak et al. [10] demonstrated a macroscopical population of ground state (in- 
CTN , plane wave number fcii) of a 2D gas of exciton-polaritons in semiconductor nanostructures (Cd-Te) at 5 K. The 
superfluid properties and Josephson dynamics of such polaritons were studied by Alodjants et al. [11] and observed 
, experimentally by Lai et al. [12]. In addition, as shown by Alodjants et al. [13] certain conditions enable optical 
cloning and quantum memory based on the cavity polaritons in question. 

In this paper, we discuss models of polaritonic crystal (PolC), which can be produced using existing technologies 
and procedures for laser control of atoms. A remarkable feature of such structures is the possibility of complete 
polariton localisation, an analogue of light localisation in photonic crystals in nonlinear optics. This effect markedly 
reduces the group velocity of an optical wave packet propagating through the medium. At the same time, polaritonic 
crystals can be used to observe BECs of low-branch polaritons. 



II. MODELS OF POLARITONIC CRYSTALS: BASIC EQUATIONS 



Consider two models of PolC. In model 1, an ensemble of ultracold two- level atoms is confined in a deep optical 
lattice (Fig. la). This can be done by a number of experimental means, in particular by using a two-component 
(spinor) condensate of atoms with internal levels \a) and \b) [14, 15]. In this case, a 2D periodic structure of elliptical 
(needle-like) atomic condensates can be produced using interference of two standing waves (not shown in Fig. la) 
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FIG. 1: Schematic representation of the two models for polaritonic crystals: (a) trapped ensembles of ultracold atoms in a 
cavity interact with a quantised light field, whose distribution is represented by grey circles; (b) polaritonic crystal formed by 
an array of microcavities containing macroscopic ensembles of two-level atoms. 

along the x and y axes, respectively. The atomic ensembles will then interact with quantised light field in the cavity 
along the z axis in the strong coupling regime (see Eq. (2) below) . 

Model 2 considers a lattice of M tunnelling-coupled microcavities in the xy plane (Fig. lb). Each cavity contains 
two- level atoms interacting with the electromagnetic field along the z axis. The polaritonic crystals differ fundamen- 
tally from coupled- resonator optical waveguides (see, e.g., [6]) in that they offer the possibility of photon tunnelling 
in the xy plane, normal to the main axis of the cavities. 

It can be shown that, in the limit of so-called strong coupling between neighboring cells (cavities in model 2) 
containing atoms, the two models schematically illustrated in Fig. 1 give the same physical results. In this approxi- 
mation, an atomic system can be described as a system of bosons (bosonic modes) that evolve only in time: its spatial 
degrees of freedom are frozen. This approximation is valid when the number of atoms N in each cell is relatively small 
(N < 10 4 ) [14]. The height of the potential barrier between the cells of the optical lattice far exceeds the chemical 
potential of each atomic ensemble [16]. Otherwise, it is necessary to take into account the spatial configuration of the 
atomic system, which leads to the formation of spatially localised atomic structures [17]. 

Here, we restrict ourselves to the strong coupling approximation, completely neglecting interatomic interactions 
and taking the atomic ensembles in the cells to be an ideal gas. It will become clear from the analysis below that the 
polariton model for atom-photon coupling is then quite correct. 

To be more specific, we consider model 1 of polaritonic crystal. The strong atom-photon coupling condition is 
thought to be fulfilled: the coupling parameter k in each cell of the lattice is substantially greater than the inverse of 
the coherence time, r co ^, of the atom-photon system: 
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k > 



(1) 



The total Hamiltonian of the system is 



H — H at + Hi n t + H p h, 



(2) 



where 




(3) 




(4) 
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Here, H at represents the ensemble of noninteracting two-level atoms in the trap; H int represents the atom photon 
interaction in the cavity in the rotating wave approximation; H p h represents the light field in the cavity in the paraxial 
approximation; are the boson annihilation (creation) operators for the levels j — a and b; M at is the mass 

of a free atom; A is the Laplace operator; V^H is the total atom trapping potential, which comprises the harmonic 
potential of the magneto-optical trap and the optical lattice potential along the x and y axes [18]; \E'(4' + ) is the 
annihilation (creation) operator for a field propagating along the z axis of the cavity, M p h = hk z /c is the photon 
effective mass in the cavity; k z is the z-axis projection of the optical field wave vector; V p h is the photon trapping 
potential in the atom-photon coupling region, which can be created by special gradient-index lenses or fibres [8] ; and 
Aj_ = d 2 /dx 2 + d 2 /dy 2 . 

For an array of M cells in an optical lattice, the $j and '5 operators can be represented in the form 

M 

$a = Yl a m(t)^ m (f), (6) 
rn — 1 

M 

MtRW. (7) 

m— 1 
M 

*=£ V>m(t)6n(r), (8) 
m— 1 

where f^{r) and £ m (r) are the real- valued Wannier functions describing the spatial distributions of the atoms and 
field, respectively, in the mth cell (m = 1, M). In the limit of strong coupling between neighbouring cells (so-called 
tight-binding approximation), the <f^ b functions satisfy the relations [16, 17] 

j \v a m "\ 2 d 3 r= 1, J tftf +1 rf 3 r- 0. (9) 

Analogous relations are valid for the ^ m (r) functions. The a m (t) and b m (t) operators characterise the dynamic 
behaviour of the two components (two modes) of the atomic ensemble at the lower and upper levels, respectively, and 
the tl) m {t) operator describes the time evolution of the cavity field in the mth cell of the lattice. 
Substituting (6)-(8) into (3)-(5) we obtain 

H at = hY^ [^mat a m a m + ^mat b m b m - y (a+a m _i + a+a m+ i + h.c.) - — (6+6 m _i + h.c.)}, (10) 

m— 1 

M 

H^t = ft 9 m {^a+ l b m + b+a m ip m ), (11) 



M 

H ph = hY^ [^mph^ti^m ~ -(^+Vm-1 + V^m+l + h.C.)], ( 12 ) 
m— 1 

where the coupling coefficients [3 a ^ and a characterise atom (photon) tunnelling between neighbouring cells and 
are determined by the overlap integrals of the (p'j n (r) and £ m (r) functions with their derivatives, respectively. The 
quantities ui^at an d &mph are defined in an analogous way [16, 17]. We take all the atom-photon coupling coefficients 
to be the same in all the cells: g = gi = g 2 = ... = gu- 

Let us turn to the momentum (fc-space) representation. Given that polaritonic crystals have a periodic structure, 
the ip m , a m and b m operators can be represented in the form 

a m = --L= £ aj:exp(ikn), b m = — L= ^ b^exp(ikn) 7 ip m = — L= ^ tpj:exp(ikn), (13) 

k k k 
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where ft is a lattice vector. 

For simplicity, we examine a ID polaritonic crystal structure below, for which kft = mk x l, where k x is the x- axis 
projection of the optical field wave vector and I is the lattice constant. Substituting (13) into (10)-(12), we obtain the 
/c-space Hamiltonian 

H = h ^ (^(fc)V^ + \" at (k) (bib, - + 0l ) + -j= £ (^ a Put + £^f)) • ( 14 ) 

Here u) p h(k) and uj a t{k) determine the dispersion relations for the photonic and atomic systems of the polaritonic 
crystal, respectively, and are given by 

uj ph {k) = LO m ph - 2acos(kl), u at {k) = cj h mat - ^ nat - 2/3cos(kl), (15) 
where (3 — fli, — [3 a is the effective coupling coefficient of the atomic lattice. 



III. QUANTUM DEGENERACY OF A ID POLARITON GAS 



In the strong coupling regime, expression (14) is a many particle Hamiltonian in the momentum representation, 
which describes a ID periodic structure and can be analysed in terms of dark- and bright-state polaritons. We will 
consider it in the low atomic excitation density limit, where all the atoms predominantly occupy the lower level \a) 
[11]. The boson annihilation (<f>~) and creation (<pi) operators for collective excitations in a two- level atomic system 

can then be defined in the fc-space representation: 



4 = E 



g i 
Using (16), Hamiltonian (14) of the system can be represented in a more convenient form: 



(16) 



(17) 



where we denote again g instead of gVN . Hamiltonian (17) can be diagonalised using the Bogoliubov transformations 

s a = - V2<j>w E 2,u = vih + M2^, (is) 

where 



2 1 

Mr, 2 = 2 



it 



6u> 



{5lj 2 + Ag 2 ) 1 / 2 



(19) 



are Hopfield coefficients satisfying the normalisation condition = ^] Suj = bj at {k)—Lo p h{k) = A— 2(/3— a)xcos(fc^) 

is the frequency detuning, dependent on quasimomentum k; and A = u h mat — cu^at ~ u mp h is the detuning at 
kl = 7r/2 + irp with p = 0, 1, .... 

The E^g and S 2 £ operators in (18) represent two types of elementary excitations in an atomic system upper and low 
branch polaritons, with characteristic frequencies li2 (fc), that determine the dispersion relations and band structure 
of PolC. The frequencies are given by the expression 



fii, 2 (fc) 



U a t(k) + LU ph {k) ± ^/SiU 2 + 4 5 2 J , 

Using (20), one can find the mass of the upper (subscript 1) and low (subscript 2) branch polaritons: 

2m at m ph A 2 + 4g 2 



(20) 



TOi2 = 



(21) 



(m at + m ph )J A 2 + Ag 2 =F (m ot - m ph )A 
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FIG. 2: Polariton dispersion branches fii (dashed curve) and 0,2 (solid curve) as functions of reduced quasi-momentum (Bloch 
vector) in the first Brillouin zone. The characteristic frequency of the atomic transition is (w^at ~ ^mat)/^ 71 = 500THz, A = 0, 
the atom-photon coupling frequency is g(2iv)^ 1 — 2GHz, the photon (atom) effective mass in the lattice is m p h — 5 x 10~ 36 fc<7 
{m at = 38.5 x 10~ 27 kg), and the lattice constant is I = 2.24pim. The inset shows the frequency range of Rabi splitting. 



where A = A — 2(/3 — a) is an effective detuning, which includes the characteristic frequencies a and /3, m at — h/ (2f3l 2 ) 
and m p h = h/(2al 2 ) are the effective masses of the atoms and photons in the lattice, respectively. 

Figure 2 shows the upper and low branch polariton dispersion, £!i, 2 (fc), in the first Brillouin zone. The band gap 
is determined by the Rabi splitting: {5ui 2 + tig 2 ) 1 / 2 . In the center of the zone (near k — 0), both dispersion branches 
are parabolic. At resonance (Suj = 0), the splitting is governed by the atom-photon coupling coefficient 2g (see also 
the inset in Fig. 2). 

The minimum at k = in the lower polariton branch in Fig. 2 is of fundamental importance. At small magnitudes 
of the quasi-momentum, kl <§C 1, we have from (20) 

hk 2 

n 2 (k) « — , (22) 

which corresponds to the dispersion law of free particles (polaritons) at the minimum in 0.2(h) in Fig. 2. The 
statistical properties of the polariton gas are then governed by its dimensionality (see, e.g., [19]). In particular, in the 
case of resonance coupling the low branch polariton mass can be found from (21): 

^ _ 2m ph 
2 l + mph/niat' 

For m p h/m a t 1 (a 3> (3), the polariton mass is sufficiently small. For example, in the case of interaction of 
two-level sodium atoms with an electromagnetic field in a cavity at a level separation wavelength of 589 nm, the 
photon effective mass in the cavity is m p h — 0.5 x 10~ 35 fcg [12]. The polariton mass estimated by Eq. (23) is then 
77i2 — 10 kg. Therefore, the quantum degeneracy temperature of a ID polariton gas, Td = 2Trh 2 n 2 /(m,2kB), may 
be rather high (~ 300-ftT at a polariton density n\ ~ 10 4 cto _1 ). It should, however, be kept in mind that the spatially 
periodic structure of the polaritonic crystal in Fig. 1 has a long coherence time only at sufficiently low temperatures, 
where a relatively large number of atoms can be trapped. The polariton gas can be considered a highly degenerate 
quantum system, meeting the condition tiiKt 3> 1, where At — [2ttH 2 /(m^ksT)} 1 / 2 is the de Broglie wavelength. In 
this limit, the formation of coherent polaritons in a polaritonic crystal is of interest for spatially distributed recording 
and storage of quantum optical information [20, 21]. 
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IV. GROUP VELOCITY OF POLARITONS 



Consider the group velocities vi, 2 — dCl\^{k)/dk of polaritons in a lattice. From (20) we obtain 



h sin (kl) 

Vl.2 = 



2lm 



ph 



m ph f m ph \ 8uj 



T 1 



m at ' V m at J yjfoj 2 + Ag 



(24) 



It follows from (24) that the group velocities t>i.2 of polaritons are low at small k. In particular, in the kl -C 1 limit the 
low branch polariton has a linear V2(k) : V2 — hk/rri2- At the same time, v\ t 2 = at the boundaries of the Brillouin 
zone, i.e., at kl — pTr 7 p = 0, ±1, .... In this case, the structure of the polaritonic crystal allows polaritons to be fully 
localised within this zone. 

The ability to reduce the group velocity of polaritons can be used to observe 'slow' light, which is of high current 
interest for quantum optical information recording and storage. In our case, the group velocity of the optical field 
can be varied by changing the atom-field detuning A (or A). 

The possibility of controlling low branch polaritons with small magnitudes of the quasi-momentum, k, is exemplified 
in Fig. 3. Control is performed near the bottom of the well in Fig. 2, where the parabolic dispersion law (22) is valid. 
In this limit, polaritons can be assigned a wave function, ^(x, t), which carries quantum information and satisfies free 
particle Schrodingcr equation: 

*4 + ^l£>W,*) = 0. (25) 



dt 2ui2 dx 2 

The solution to Eq. (25) is well known in quantum physics (see, e.g., [13, 22]). Polaritons are a coherent wave 
packet that broadens with time and propagates through the polaritonic crystal (Fig. lb). The characteristic wave 
packet broadening time, t = 1112 f 2 /h, depends on both the polariton mass, TO2, and the x-axis width of the packet, 
/, at the initial instant, that is, on the incident beam diameter. 

Quantum optical information recording and storage using PolC is expected to have a three-step physical algorithm 
based on the ability to control the group velocity of a polariton wave packet in the medium by varying the A detuning 
(Fig. 3a). 

For quantum information to be recorded in step I, the A detuning must fulfil the condition A 3> 2\g\. The 
corresponding time interval in Fig. 3a is < r < 0.25. The low branch polariton is here photon-like, i.e., S 2 £ ~ — tp^ 
(fii ~ 0,yU 2 — — 1), and its effective mass m 2 — m p h [see Eqs. (18), (19) and (21)]. The wave packet then propagates 
at a velocity given by 

v 2 = — = 2al 2 k, (26) 
m ph 

and displaces along the x axis from the point x = to x = v 2 t (Fig. 3b). For example, at a quasi-momentum 
k = 10 5 m _1 the estimated polariton group velocity is V2 = 2 x 10 6 m/s. 

To map the optical information to coherent excitations of the medium in step II, the detuning should be made 
negative: A <C — 2\g\. In this limit, the low branch polariton becomes atom-like, so that S 2 £ ~ (/ii ~ l,/i 2 — 0) 
(Fig. 3). At m p h/m at — j3/a -c g 2 /A 2 , its group velocity is [7] 

hkg 2 2al 2 kg 2 

v 2 = ^- = — (27) 

m ph A 2 A 2 

When the more stringent condition -c <C 1 is fulfilled, the low-branch polariton group velocity is 

v 2 = — = 2(3l 2 k, (28) 
m at 

which corresponds to the velocity of the atoms in the lattice. In particular, for sodium atoms with an effective mass 
m a t — 38.5 x 10 _27 fc.g the group velocity of such polaritons at the above wave vector is V2 — 2.6 x 10~ 4 to/s [21]. 

In fact, Eq. (28) determines the lower limit of the velocity of an optical wave packet in the structure of a polaritonic 
crystal at a given quasi-momentum. In this limit, all the information carried by a light beam is recorded and stored 
via atomic excitations. At an incident beam diameter / ~ 10~ 4 m, the estimated characteristic broadening time is 
T5 = m a tf 2 /fo ~ 3.7s. This time scale (the interval 0.25 < r < 0.75 in Fig. 3a) determines the longest information 
storage time in a quantum gas of two-level sodium atoms. More precisely, a necessary condition for such information 
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FIG. 3: Propagation of a ID polariton wave packet through a polaritonic crystal: (a) D — A/(2\g\) as a function of normalised 
time r = ht/(m p hf 2 ) and (b) wave packet envelope (probability density) S = |\E'(a;,t)| 2 /|\E'(0, 0)| 2 as a function of r and 
normalised coordinate X = x/f at k x f — 10; optical information (I) recording, (II) storage and (III) reading (retrieving) steps. 



recording is t st0 r -C n (where t stor is the information storage time in the atomic system), which implies that the wave 
packet retains its shape during the whole period of quantum state storage. 

To restore (read) optical information at the output of the medium after time t s tor (step III), the polaritons must 
be again made photon-like by switching the A detuning in the reverse direction. In Fig. 3a, the characteristic reverse 
switching time of the wave packet, t retr , is determined by the time interval 0.75 < t < 0.8. At the output (r = 1), 
we again have an optical wave packet, which is displaced along the x axis in a plane normal to the cavity axis. 

In this work, we do not assess the quality (fidelity) of information storage. Alodjants et al. [13] estimated the 
information storage fidelity taking into account only changes in wave packet shape, but this is generally insufficient. 
It is, in addition, necessary to analyse transformations of the quantum state of the optical field with consideration for 
the structure of the polaritonic crystal and its decoherence (see, e.g., [23]). Analysis of this problem is of interest on 
its own and is beyond the scope of this paper. We note only that, under real experimental conditions, the information 
recording, storage and reading time is limited by the decoherence time of the polaritonic crystal. It is therefore quite 
reasonable to use a condensate of atoms that have a sufficiently long macroscopic coherence time-tens of microseconds 
according to experimental data [24] . 



V. CONCLUSIONS 



We examined a lattice model of coherent polaritons in a spatially periodic structure a polaritonic crystal formed by 
a lattice of ensembles of two-level atoms effectively interacting with quantised electromagnetic field in a cavity (or in 
a 2D lattice of cavities) in the strong coupling regime. Our results demonstrate that the structure of PolC allows the 
low branch polariton to be fully localised, which can be used, first, to achieve quantum degeneracy of the polariton 
gas and, second, to substantially slow down group velocity of light pulses in such media. The coherence properties of 
an ensemble of polaritons were discussed from the viewpoint of spatially distributed quantum recording, storage and 
retrieving of information related to a propagating optical wave packet. 
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